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ABSTRACT 



This investigation is a study of the problem faced by a circuit 
designer who In the analysis and synthesis of his circuit must use 
parameters which vary over wide ranges « 

A technique utilizing a digital computer resulting in design 
curves is developed. A second technique is demonstrated where the 
same design curves are obtained by another method employing an analogue 
computer. 

Both techniques are applied to a neon pi ot oconductor circuit and 
design curves are obtained with each method. 

The writers wish to express their appreciation for the assistance 
and encouragement given them by Professor W. Conley Smith of the U. S. 
Naval Postgraduate School, 
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1. list |■©d^sc£ loti 



Sysfigi® reliabtlifcy is th@ keynote in this >^yin©'mle ag@ of saeehaniz- 
atlon wh®ra unfailing performance is expected and required. However, 
component parts in a systenas configuration are often subject to manufactur- 
ing tolerances which raay very radically from the norm, yet a systems de- 
signer has to deal with these varying components and produce a dependable 
system. Dependability can be measured as a statistical quantity or as an 
absolute requirement. The latter definition of dependability will be used 
in this paper as a yardstick In measuring the success of a particular de- 
sign, The problem to be undertaken in this thesis is the synthesis of a 
reliable system consisting of widely varying components. Reliability as 
used here is dependent on the components themselves not falling. 

Circuit synthesis is generally Interpreted as the formulation of sig- 
nal flow charts and the determination of component values for the realiza- 
tion of a given transfer or immlttance function. This paper will be con- 
cerned with only the latter part of this definition. 

The specific system to be investigated and designed consists of a 
neon photoeonductor logic circuit which is a basic building block for other 
logic circuitry. Justification of the Important requirement of reliability 
in this circuit is indicated by Its use as a basic element In the founda- 
tion of an entire system. Analysis and synthesis of the circuit will be 
restricted to a steady state investigation. 

Using the specific neon photoconductor logic circuit as an example, 
a general technique of circuit synthesis will be proposed in which the 
tools of circuit analysis and Inductive reasoning are utilized for synthesiz- 
ing. Analysis can exist without synthesis but the converse of this state- 
ment is not true, H. S. Scheffler and F. R, Terry in their paper entitled 
"Description and Comparison of Five Computer Methods of Circuit Analysis" 
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descrife® three non-statistteai approaches t© circuit analysis, Sonie of 
their concepts in circuit analysis are adopted and extended Into the field 
of circuit synthesis. For example the criteria, point of failure, is a 
basic extreme condition applicable to a circuit problem whether the investi- 
gation be of analysis or snythesis. The particular logic circuit to be 
investigated will be synthesized with the aid of both a digital and an 
electronic analogue computer. Each computer has advantages over the other 
which are peculiar to the methods utilized in the solution. Successful 
circuit synthesis by the computers produces design curves, which when 
properly Interpreted can be an Invaluable asset to the circuit designer. 
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2« Descript JI. 0 R of Digital C<»p«Jt©r Design Technique® 

Ihif s^scfciaa ©xplalssS;, tn general CeruiSj s technique useful in cir- 
cuit design which uses to full advantage the high speed calculating feature 
of a digital cemputero 

In designing a circuit to tn®et output specifications the designer is 
alrnost invariably faced with the problem of using eomponents In the cir- 
cuit which do not have constant values*, but may have values in a given 
range. The simplest example of this point is an ordinary fixed resistance 
which ffiiay have any value in a range abov'e and below its nominal value. 

This variation is due to manufacturing differences, A range of variation 
may also be duo to environmental changes during circuit operation^ aging 
or other factors. 

Circuit parameters which must be selected from within a known range 
of variation will be referred to as "'known" parameters. Other parameters^ 
whose value can be chosen by the designer based on his design to meet 
specifications are defined as "design" parameters. 

The final circuit or circuits decided upon by the designer must meet 
the specifications for any random set of the "known" parameters - each 
parameter in the set picked arbitrarily from anywhere in its range of 
variation. The designer must choose the design parameters for Ms circuit 
so that when "acceptable" values of the design parameters ar® used In the 
circuit, the circuit performance will meet specifications for any random 
set of the known parameters. 

In simple circuits the decision as to accepteble values of the design 
paramters may be easy. Consider a simple volcage source and resistance 
circuit. If the voltage is a known parameter with variation from 5 to 10 
volts tmd tha specification on output demands a current of fror. 3 to 4 
amperes, ihet» the selection of acceptable values of reslstanc®j, the design 
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parnu'u'‘jt er i ^ s * ^<)||.> ■ , 

A . - * 
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If the design parameter^ R, hc S any aceept^^ble yalu^; that Is am- 
value from 1.25 to 3.33 ohms, the circuit will meet perfomance specifica** 
tions for any value of the knovn pAramcter, '>oltage, Ir^ its range 5 to 
10 volts. 

The example shows the significance of the terras known and design 
parameters, but is much too simple to demonstrate the technique developed. 
In the general case, there are many knowi parartjeters, each with a range of 
variation. There is usually a number of output specifications to be met 
and the designer is at liberty to pick his design parameters in sets of 
acceptable values. 

Stated most simply, the technique begins by first writing equations 
which express the output as a function of the design and kno^^ parameters. 
By inspecting these equations and employing boundary conditions, the de- 
signer must establish a range of Investigation for each of the design 
parameters. This range is similar to the range of variation of the knowiji 
parameters, but differs in that it may be made larger or smaller during 
the process of design based upon the designers interpretation of the in- 
itial answers. 

Once a range has been established for each parameter, the equation 
expressing the output must be solved for each value in the range of each 
parameter. The values of design parameters which yield acceptable output 
regardless of the values of the known parameters used are recorded and 






plotE® '‘■r <]esl8,R pt'Cflraettir^™ a c4te«ie 

feutlt wfth an aecapfci^bld earabte^ition ©f chess desigr- j.-ar9»>eE®?? will pro- 
vide oytpMt which meets specif Icatloiag no matter vhist valwes of knoim 
paranvietars are used* 



Since It is Iraposslble to solve the circuit equation for averv value 
of every design parameter end every known pararaeter, the ranges of th«isa 
parameters must be divided Into increments and the equation solved for 
each Increment @f each parametero The selection of theso increments in- 
troduces sor^ error* As in other approximation procedures, smaller in- 
crements introduce less error, Ho^^ever, smaller increments lead to is-cre 
coKibtnatlons to investigate and thus many more celeulations. If is 
the number of Increments of the n th parameter t© be investigated, the 
total number of answers to b@ computed and tested Is 



It is seen that any reasonably complicated design problem with a 
practical number of IncreiRgnts quickly leads to an extremely large sus^feer 
C'f cstculacicia.^ , This problem suggests use of a digital computer with 
its rapid calculating ability as an invaluable tool. 

In some cases it may not be necessary to investigate all increments 
of a known parameter. This exception was demonstrated in the trivial ex- 
ample seated above although not specifically pointed out. The range of 
acceptable resistance was determined by investigating only the maximum 
and minimum values of the known parameter, voltage. Because of the simple 
nature of the circuit It is obvious that a value of the design parameter 
which is acceptable at both the maximum and minimum values of the known 
parameter will be acceptable for any value between and only the end points 
need be used. The characteristic of the parameter voltage which makes 





<•; J, 



is always of the same sign. 



this true is that the partial derivative 
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Trs Cher*, it the partial derivative of the output taken with 

respect to a known parameter is always of the same sign» for all values 
of all the other known and design parameters^ only the maximum and mini” 
mum values of that known parameter need be Investigated^ since if the 
design parameters are selected to be acceptable with both these points, 
they will be acceptable for any value of this parameter In its range. 

Since parameters of this sort mean a great reduction in the number 
of calculations, the equations for output should be partially different- 
iated with respect to each known parameter in turn to determine if deriva- 
tives of constant sign exist. The partial derivative must be of the same 
sign for al 1 values of the other known parameters. This may sometimes be 
shown by manipulation of the derivative expression to show It is always 
of one sign due to Its form. In other cases a worst case technique will 
be sufficient to prove a constant sign of the derivative. In this pro- 
cedure, if the derivative is asstimed to be always positive all known para- 
meters appearing in the expression which tend to increase the derivative 
are set at their minimum and all those which tend to decrease the deriva- 
tive are set at their maximum. If the derivative remains positive for this 
worst case. It can be assumed always positive. This procedure has the 
weakness that It may not be evident whether a given parameter tends to 
Increase cr decrease the derivative. Also, in general. If the range of 
parameter variation is very large, the derivatives will not pass this 
worst case test and the parameter must be divided into Increments for in- 
vestigation. 

When the ranges of all parameters have been divided Into increments 
to be tested whether numerous Increments or only one (end points type of 
parameter) the computer program Is written. This program solves for the 
output using each increment of each design parameter in combination with 
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aach isi- r^jnertt o£ each known pataineeer. The output resulting from each 
set of parameters is tested and those sets which result in output which 
meets specifications are stored and printed. 

The output of the computer can most easily be interpreted if pre- 
sented in graphical form. If there is Just one design parameter this 
is not necessary, A most useful presentation for two design parameters 
is a plot of one versus the other. Starting with the data for one value 
of a selected known parameter; for each value of one design parameter, the 
maximum and minimum acceptable values of the other are plotted. Connect- 
ing all such points will yield an area. The values of design parameters 
corresponding to any point within this area when used with the previously 
selected known parameter will result in a circuit which will meet design 
specifications. 

This procedure is repeated until area plots are obtained for every 
Increment of each of the known parameters. In each case one design para- 
meter Is plotted versus the other utilizing the same scales. If all these 
plots are overlayed, the area which is common to the acceptable area of 
all the graphs is the final design area desired. If a point In this final 
common area is selected and the corresponding values of design parameters 
are used in a circuit, the circuit will meet specifications for all known 
parameters. 

The requirement for obtaining an area graph for each increment of 
each known parameter demands much tedious plotting. This is essentially 
the procedure; however, it is possible to write the computer program so 
that the data will be partially processed by the computer before printing 
thus eliminating some of the plotting and overlaying. This programming 
method Is most easily described by use of a specific example and is de- 
monstrated in a later chapter. 
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Ir '.B quRts. posstfc'iw that when the acceptable area graphs for the 
incretnents of, known paran'ieters are overiayed there will be no acceptable 
area which is common to all the graphs. Thus there are no combinations 
of design parameters which will guarantee satisfactory circuits for ran- 
dom known parameters. This condition can occur because the specifications 
on the output are too strict. This may be stated another way by saying 
that the range of one or more of the known parameters is too large. The 
difficulty may be corrected by relaxing the specifications on the output 
or dividing one of the known parameters Into two groups. If the trouble- 
some known parameter is divided Into two groups, each of which has half 
the range of variation of the parameter, there may then be common areas 
of acceptable values for these half ranges and all other known parameters. 
In this case then there are in effect two design areas rather than one and 
the divided parameter must be tested and separated into two groups before 
one can use this procedure to design a circuit. 

If the area is large, it may be possible to tighten the output speci- 
fications to be met or to add another specification. Doing either of 
these will tend to reduce the size of the area. By inspecting the size 
of the area resulting from the initial output specifications, a decision 
can be made as to how the specifications can be altered. By Inspecting 
subsequent results and making further adjustments a balance between specif- 
ications and available combinations of design parameters can be achieved. 

The technique described results in an area of acceptable design 
values. If the designer picks a set of design values from in the final 
area, he is assured that the circuit i#!! t maet oatput sp^joificatlons. 

The question still remains, "If all the combinations in the area meet 
specifications is there not one best combination or some part of the area 
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^vhirh circuits than otafer parts?^^ The tecl*miqy@s des- 

cribed up to the point will not answer this question, 

t)ne ftiethod of determining the best part of the final common area Is 
to select one of tha preliniinary area plots. This plot corresponds to 
one set of knora parameters. Using this set of known parameters and re- 
presentative values of design parameter combinations throughout the area, 
one calculates the value of one of the output quantities. The values cat 
culated are plotted at the corresponding point in the area. Inspection of 
the plot may reveal a trend in the value of the output parameter and In- 
dicate a portion of the area where this value best approximates Its opti- 
mum value. This process may be repeated with other sets of known para- 
meters and the corresponding area until It Is established that the trend 
holds in enough of these prelirriinary areas that It may be assumed to hold 
in the final area plot. 

In sumaitary tha technique which has been described in this section 
results in a plot of one design parameter versus the other. An area of 
acceptable combinations of the two design parameters can be drawn on this 
plot. If a combination is selected from inside this area and a circuit 
constructed using this combination and any combination of the known para- 
meters used in the circuity the resulting circuit will meet the desired 
output spec I f lea t ions . 
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3* Desctlfitlon of ^v.o^-*photcccndlJlCtot Circuit* 

The circuit chosen to illustrate the technique described briefly In 
the preceding chapter is a neon tube * photoconductor circuit which is 
used as a basic building block in the design of computer logic circuits. 
The basic circuit is shown tn figure 1* 



r 

) 




Figure 1, 

The circuit consists of a neon tube, two fixed resistances^ R1 and 
R2 and a photoconductor^ Rpl, supplied from a 300 volt source. 

In the circuit, the output quantity of interest is the light of neon 
tube two. The input is light from a similar external neon tube, called 
number one. With no Input, there is no output; the very high dark resist- 
ance of the photoconductor which is in the mega ohm region limits current 
flow and keeps the voltage drop across neon number two low enough to pre- 
vent its firing. With input light from neon number one,^ the resistance 
of the photoconductor drops radically to its light resistance value of a 



few kllohms amid neon oufrtber two fires^ producing light output* 

The light output from neon number two is used as input to sub* 
se circuits so it is essential that it be as nearly constant as 

possible for stable operation* This constant current characteristic 
is not possible with the circuit as shown since the resistance of the 
photoconductor is a function of temperature. As the circuit operates j, 
the ambient temperature riseSj causing the photoconductor resistance to 
rise and the neon current to decrease* A proposed method of compensation 
intended to minimize this temperature effect is to insert a similar photo 
conductor in the shunt branch of the circuit as shown in figure 2* 
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Figure 2. 

In this compensated circuit the output current is more nearly con- 
stant despite temperature Increase because of the feedback nature of the 
output light shining on photoconductor Rp2* Now^ as the output current 
tends to decrease as a result of temperature variation changing the re- 
sistance of Rpl, the light output decreases^ Increasing the resistance of 
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Rp2 and shtmt ing ot the aval ible current through the neon tending 

to maintain the current constant. 



The problem then is to design a circuit, using available components, 
which will give a nearly constant output light or current, throughout 
the operating temperature range of 25 to 55 degrees Centigrade. The 
supply voltage is assumed to be a constant 300 volts. The neon to be 
used is a General Electric NE2 type which has a firing voltage of 135 volts 
and a maintaining voltage of from 60 to 100 volts. 

The characteristics of the photoconductors employed in the circuit 
are more complex since the value of photoconductor resistance is a func- 
tion of both the ambient temperature and the incident light. The latter 
variation is more easily investigated if neon current rather than light 
is considered the determining factor. Mathematical relationships governing 
photoconductor resistance variation with temperature and with current have 
been determined from experltrental data. 

The temperature dependence of resistance can be obtained by consider- 
ing figure 3 where 1°8 jq resistance in kilohns is plotted versus tempera- 
ture in degrees Centigrade for three typical photo conductors. The curves 
are all for a constant incident light corresponding to two milliantperes of 
neon current. Since the data plots as a straight line on semi- log paper 
the variation Is of the form Kp" /n c. , if curves II is selected to 
determine the constant. 




f K / C ' A' 



( 1 ) 





( 2 ) 



Dividing (2) by (1) 
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Thus 









12 



1 

I 



1 



f/iLoHiHS 




3 





i 







I 



The coefficient A and the mratheirtat leal relationship expressing re- 
sistance variation \vith neon current can be obtained by considering 
figure 4, In this curve^ photoconductor resistance in kilohms is plot- 
ted versus current in the exciting neon expressed in milliamperes ; the 
curves are all for a constant temperature of 40"^ C. The thr^^e curves 
kf B, and C are for t^ree representative photoconductors ^ A and C show 
the extremes of variation and curve B is considered typical. In attempt- 
ing to express the photoconductor resistance as a function of light, it 
was noted that curves A, B and C are all of the same mathematical form 
and are related by a constant. This nmgnitude multiplying factor will 
be defined as Z. This relationship is best illustrated by table I, figure 
5. 

The current-resistance characteristic of any photoconductor can be 
obtained by deriving a mathematical expression which will fit one of the 
curves giving the proper form of variation and then multiplying the en- 
tire expression by a constant to fix its magnitude. 

Curve B, selected as typical will be the one used to determine the 
mathematical form of variation. If this expression is assumed to have a 
magnitude multiplying factor of unity, reference to figure 4 shows that 
in order to cover the range from curve A to curve C the magnitude multiply- 
ing factor must have values of from ,4, for curve A to 2.5 for curve C . 

The first approximation to curve B which was investigated is a 
hyperbolic function where is a constant to be determined. 

^ r I 

Since a neon current of 2 ma. will be used as the mean point about which 
to design, the hyperbola constant was determined to fit curve B exactly 

3 

at that point. From tie curve at I ^ 2.0 ma., Rp ^ 4.388 x 10 ohms. 
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Table 1 



Current 
in ma. 


Rp 

Cuir/e A 
in Kohm 


Rp 

Curve B 
in Kohm 


Rp 

Curve C 
in Kohm 


Rp(A) 
Rp (B) 


Rp(C) 
Rp (B) 


1.00 


4.48 


12,66 


30.8 


.354 


2.43 


1.25 


3.38 


9,07 


22.0 


.370 


2.43 


1.50 


2.78 


6.71 


17.0 


.414 


2.53 


1.75 


2.06 


5.32 


13.86 


.388 


2.60 


2.00 


1.72 


4.388 


11.22 


.394 


2.54 


2.25 


1.46 


3.79 


9.28 


.396 


2.44 


2.50 


1.26 


3.37 


7.83 


.374 


2.34 


2.75 


1.11 


3.04 


6.76 


.363 


2,22 


3.00 


1.00 


2.89 


5.82 


.346 


2.02 


ki 


- = 1 


Figure 5 
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Incorporating the effect of temperature, 

i: 7 




3 

Since the value of R = 4«388 x 10 ohms from curve B corresponds to 



I * 2ma. and T si 40 ®C the constant 8*776 K 2 can be determined 






J X } 

f A-;2 ^ -^Z- ^ ^ 

The hyperbolic expression with temperature effect included approxlmat ing 



curve B is then 



)?,. 



- - c 

X 

To include the fact that the magnitude of the resistance can vary 
from curve A values to curv€ C values, the magnitude multiplying factor 
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Z mast be included. Tho complett; resistance expression for the photo- 
conductor Is then 



In the circuit being analyzed, Rpl is supplied with exciting light 
from an external neon, neon number 1, which is assumed to provide a con- 
stant light intensity corresponding to a neon current of 2 ma. There- 
fore; I " - ' 




/ 






/ 



c / -f ^ 



The compensating photo conductor Rp2 Is excited by light fed back 
from the output neon, neon number 2. This light corresponds to the out- 
put current Therefore: 

2 _ . c y ^ 

aV ^ 

z ^ 

The plot of the hyperbolic approximation to curve B shown in Figure 
6 indicates that the hyperbolic curve fits experimental data exactly 
only at I w 2,0 ma. Since circuit performance with currents from 1,5 
ma, to 2.5 ma, is of interest, the use of the hyperbolic approximation 
introduces some error especially at the lowest value of 1,5 ma, where 
the error is 12,8%, 

The expression for photoconductor resistance is used in subsequent 
calculations and the simplicity of the hyperbolic expression Is an ad* 
vantage when used in more complex expressions. However^ a more accurate 
curve fit over the entire 1,5 to 2,5 ma, current range was desired. 

The second approximation to curve B which was investigated was de« 
rived by assuming the curve to be of the form: 

.A I RX 0 fp i~ n - ^ 

where A, B,C, and D are constants to be determined. Appendix I shows the 
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C^lcul^t ICJCilBg to ths f‘''CTl; 

. 7 ; - - V 

from which ^ 

j 1 - 

.* ■■''■ “, ■ i ' ,. i ^ 

Including the effect of temperature as in the hyperbolic approxi- 
mation: / , f ,, Y 

< f * 1 5 .J » ^ “ 

V . H i i . 

i ' ; \ “ 

In the circuit being analyzed, R1 is supplied by light from neon 
which has a constant current of 2.0 ma. therefoi^s; 



resistginc^ have such ranges: of variation, an o«Jtpat current of from 

1.5 to 2.5 will meet specifications. R1 and R2 nrnst be chosen such 
that this specification is imt w'lth any neon tube and any photoconductor 
having charactarist ics within the ranges discussed above. The only assurnp* 
tion being that both photocondoctors have the same value of magnitude multi® 
plying factor, Z. A further speci f Icat ton is that for any finished circuit 
which will have one value of neon voltage, V, and one value of Z, the vari® 
ation in current due to temperature change from 25®C to 55®C shall not be 
greater than .15 ma® The effect of Imposing stronger restrictions on the 
output will be investigated and the results compared. 

The upper limit on the values of R1 to be investigated can be estab* 
lished by setting the current II ^ 1.5 and the voltage V ® 50, 
then 

•V. - ' — 

ft . . 

R1 , camiot be as conveniently set since its value depends upon the value 
min 

of R2 which has not yet been deterunlned. A rough initial estimate can be 
obtained by assuming no current flow through the shunt branch and ignoring 
the contribution of Rpl tn comparison with R1 then 
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In view of the approximations made here, a more pessimistic value of 20 

ktlohms should probably be chosen. Experience has shown that a range of 
3 3 

from 40 x 10 to 130 x 10 ohms for R1 will adequately cover the possible 
values of R1 to have the circuit meet output specifications. 

The maximum value of R2 is indetermlnant since the effect of increas- 
ing R2 Is to shunt more current through the. neon tube. It can be seen that 
If II Is 2.5 ma. or less, R2 can be Increased without bound and the current 
I2 will never exceed the output specification. The minimum value of R2 is 
also indeterminate. In the uncompensated circuit a minimum value of R2 can 
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be deteritiined ^ lower limit exists below which the neon voltage 

will not be above the 135 volts necessary for firings In the compensated 

circuit however;, dim to the very high dark resistance of the compensation 

photoconduccor, the neon will fire for AX\y value of R2. The mcignltuda 

of R1 . was selected as an initial lower limit for R2 also* An initial 
min 

upper limit of 700 K ohms was selected* Subsequent experience showed 

3 3 

that a range of R2 from 40 x 10 to 530 x 10 ohms was sufficient to yield 
the desired design curves. 
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4 . P re 1 Iml ria Ty i rc n l lu km^ ? y s !. . 

This secctoim deals with the application of the de-sign technique deva» 
loped in Section 2 to the neon^photoconductor circuit which ha.s been pre» 



vioualy descriHad* 

i'- 




Figure 7* 

The output quantity of interest is the light from neon number 2, 

As previously noted, the output is more easily expressed as the currant 
through the aeon, an equally good measure of performance* The technique 
first requires an equation expressing the output as a. function of the 
other parameters* 
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Then by the quadratic eqtiation: 
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iQspaction of the coefficients A,i, and C shows that since A Is 
always positive and C is always negative (the oaxifsiuee value of V is 100) 
the quantity resulting from evaluating the radical will be greater than 
B« Since current I2 cannot physically be in the direction opposite to 
that assusied, a negative answer for I2 has no significance and the equa-> 



tlon to be solved with the eoaputer is 

t / ,?y 

If the second order approxiostlon for photoconductor resistance is 



used: 
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This jsqaation can also inatiipnlated into the q^aadrattc form. 
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where the coefficients are somewhat »noxe complex, 
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The same quadratic formula must then ba solved; 
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In order to determtfie whether any of the known paranet>n‘s fall in- 
to the category where only the maximum and minimum values need to be in- 
vest igated» partial derivatives of I2 with respect to V, Z, and T must 
be taken and checked to see if they are of constant sign. 
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For the hyperbolic case: 
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Inspection Chat the quantities ‘ — ■'■ and — — are always 

po.’^itrve^ C and ara nc^gativa and the sign of B be plus 

or minus depending upon the r.t.^guli'^de of iho ni«gatiV'e term ’ ■•.. ? .■ '- ^ 

Compared vlth the i.'Of “ti.ve terai 

.. /; <"■ , '• ^ -• ' X ^ i 

The large number of terms in the expression make It unvieldly 

and no nanipulatton which would establlch tha sign of ‘ - suggests 

4 1 

itself. Therefore, the second approach must be used. The slope will be 
assumed negative and the "worst case" values of PI, P.2, T, Z, and V which 
will tend to make the slope positive will be Inserted. The expression with 
these values inserted Is: 
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The expression evaluated is a + .0374. A positive answer shows that 
the initial assumption of a negative slope may be incorrect. The infoimia- 
tion gained from this test is not conclusive, since It is Improbable that 
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the worst values of ali the v;i1I occur simultaneously. 

However^ since a positive slope ims calculated^ rlie pessinlstic alter- 
riatlve niust be cboson aud the temperature raust be investigated in In- 
crements. A similar warst case analysis on leads to the same 

conclusion. 

The variation of 12 with V will he Im^'est igated next. 
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Sau-. fc. A aiiu i' rhi.? last IneqiL^llty holds and thus the ln= 

itial assun^ptUra of a negati>/e ^ is true. Since ' is always 

negative, only the ^ximuTn and minimum values of V are investigated* 

I/lien a similar analysis of these partial derivatives is made using 
the more accurate second order approximation, it is found that , 

JjL : and '' all have possible changes in sign and 

J .) ' 

therefore v/hen this approximation is used all three parameters must be 
investigated using increments. This condition is then one disadvantage 
to using the second order approximation* 

When it has been determined that a parameter must be investigated 
using increments, the magnitude of the Increments must be chosen. As 
pointed out in Section 2, the number of increments selected is a com- 
promise between accuracy and keeping the number of calculations reason- 
able. The initial choice of increments was as follows: 

T = 25, 35, 45, 55, 

Z .4, .9, 1.4, 1.9, 2.5, 

V = 60, 100, for hyperbolic approximation^ 

V = 60, 70, 80, 90, 100, for second order approximation, 

R1 « In 10 kilohm Increments, 

R2 s in 20 kilohm increments. 

Plots made using these increments showed that parts of the range of 
R1 caused the value of R2 to increase rapidly with possible inaccuracy in 
the plot. The reason for this rapid increase is explained in Section 6, 
To eliminate this possible Inaccuracy the increments of R1 were altered 
so that in the questionable region the increments were taken much smaller 
(2 Kilohm increments). 

The problem has at this point been reduced to the solution of a 
quadratic equation with coefficients made up of the five parameters Rl, 
R2, T, Z and V, The equation Is to be solved for all combinations of 
each of the increments of the five parameters. 
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5. Digital Computer Program Technique 

The solution of a quadratic equation whose coefficients are func- 
tions of five variables and where the number of solutions is Che product 
of the parameter increments is ideally suited to a digital computer* A 
parameter increment Is defined as the number of values that a variable 
will be assigned spanning its entire range of variation. The number of 
combinations of parameters to be examined in this case will be denoted by 
N where: 









I 







Programming of a digital computer in machine language requires an 
elaborate and complex sequence of Instructions to specify every step of 
an operation. With the advent of the FORTRAN processor in late 1956, 
the intricacies of programming were no longer a stumbling block. Al- 
though the FORTRAN system was developed by J. W. Backus at IBM, it has 
found universal acceptance for all types of digital computers. The FORTRAN 
processor has made it possible for scientists and engineers to apply the 
language of ordinary mathematics to their individual problems. 

Quoting from a paper prepared by Backus and his associates, *'The 
purpose of FORTRAN was to reduce by a large factor the task of preparing 
scientific problems for IBM’s next large computer, the 704. If it were 
possible for the 704 to code problems for Itself and produce as good pro- 
grams as human coders (but without the errors), it was clear that large 
benefits could be achieved. For it was known that about 2/3 of the cost 
of solving most scientific and engineering problems on large computers 
was that of problem preparation. Furthermore, more than 90 per cent of 
the elapsed time for a problem was usually devoted to planning, writing 
and debugging the program. In many cases the development of a general 
plan for solving a problem was a small job in comparison to the task of 
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devising aas (r^chiti-a proceJures to carry out die plar*, TM goal 

of the ^ORTRAN project was to enable the prograinmer to specify a num- 
erical procedure using a concise language like that of mathematics and 
obtain automatically from this specification an efficient 704 prog^iiia 
to carry out the procedure. It was expected that such a system would 
reduce the coding and debugging task to less than one” fifth of the job 
it had been," 

FORTRAN (derived from formula translation) is a language that great- 
ly resembles the language of ordinary mathematics. Symbols such as, **, 

*> /» +> “» FORTRAN are digested as exponentiation, multlpllcaf'on, 
division, addition, and subtraction. Mathematical operations such as 
square roots, logarithms, trigonometric and hyperbolic functions are 
available in the computers library of routines and used nerely by re- 
ferencing the appropriate operation by abbreviated name. 

The Control Data Corporation 1604 digital computer was married to 
the FORTRAN processor In late 1961. The results of the merger produced 
a high speed computer that could be cosnmanded with relative ease. Suc- 
cessful results with FORTRAN programming can best be achieved by follow- 
ing these basic rules In the order shovm. 

1. Block diagram - (flow chart) 

2. Mathematical and logical operations to be performed 

3. Transformation to FORTRAN language. 

An example of a basic flow chart for the solution of the quadratic 
equation using the hyperbola approximation is illustrated in figure 8. 

The idea behind this program is to see the effect on the neon current as 
the temperature changes for any one set of known and design parameters. 
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Fig. 8 
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The fiov chart indicates the sequwica of operations to be performed 
with alternate paths iepandent on logical tests* Each block in the flow 
chart may represent a few calculations or logical interrogations. Es- 
sentially the flow chart in figure 8 systematically takes every possible 
combination of circuit parameters and tests the output current to see if 
it meets specifications. If a failure occurs the design parameter is 
altered and the calculations and logical tests are started over again. 

When the entire range of temperature variation has been investigated and 
yielded an acceptable circuit, the design and known parameter set which 
was responsible for the acceptable circuit for the entire temperature varla 
tion is recorded on magnetic tape. In addition to the parameter set the 
variation in neon current is also placed on the magnetic tape. This pro- 
cess continues with a new parameter combination until all parameter combine 
at ions have been investigated. 

The mathematical calculations to be performed in the solution of a 
quadratic equation consist of evaluating the coefficients which are vari- 
ables dependent upon the parameter set under consideration and then using 
the quadratic formula to evaluate the neon currant. Since in the present 
problem the argument of the square root is always positive and greater 
than the coefficient B, only the positive root need be calculated. T-ogl- 
cal tests to be performed consist of determining acceptable circuit designs 
by employing the given circuit specifications as test criteria. Other logl 
cal tests are used in determining maximum - minimum values and program loop 
completions. 

The transformation from mathematical to FORTRAN statements is accom* 
plished next. Input to the 1604 digital computer can be by either punch- 
ed cards or magnetic tape. Utilization of punched cards Is superior at 
this time as cards which are in error can be easily replaced^ and addi- 
tions or deletions made as necessary. The co«)plete program in FORTRAN 
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appears as in figure 9o Lach line in the program represents a single 
punched IBM card. If this same program had been done in machine Ian* 
guage, these thirty odd FORTRAN cards would be replaced by about five 
hundred machine language cards. 

If the FORTRAN program when placed in memory is complete and logi- 
cal, the program will compile. If compilation is not accomplished an 
error printout will result. The error printout categori 2 :<es the errors 
that have escaped detection by the programmer and lists these errors by 
FORTRAN statement number, Ivlten the error is recognized altetaticns to 
the punched cards can easily be accomplished and the program compilation 
retried until successftil. Writing output on magnetic tape is a separate 
entity from the computer proper and occurs simultaneously with computer 
calculations. The magnetic tape In turn is read out on an off line 
printer so as to conserve computer time. An example of the digital com- 
puter printout as obtained from the off line printer appears as In figure 
10 . 

prior to proceeding with these digital computer answers a mathemati- 
cal check is required to see whether the problem is being solved correct- 
ly, A desk calculator should be used rather than a slide rule to obtain 
enough significant figures for comparison purposes. When a calculator 
check has been made for one set of parameters and dependability of the 
computer program proven, then the accuracy of the remaining output can be 
relied upon. A similar FORTRAN program for the uncompensated neon photo- 
conductor circuit is Illustrated in figure 11. This program is needed 
for a comparison of its results with those of the compensated circuit. 
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PROGRAM BASIC 

DIMENSION V( 3) ,R 1 ( 10) ,'ft2 ( 1 7) ,T(4 ) ,Z( 4 ) 

100 F0RMATI3F4.0) 

READ 1001 V( I ) . 1= 1, 3) 

200 FORMAT ( 10F7.0) ,,, , 

READ 200 IRl (J) , J= ), )6l 
300 FORMAT ( 11F7.0) 

READ 300(R2(K),K=1, 17) 

400 FORMAT (4F3.0) 

READ 400 (TIL), L=1,4)' 

500 FORMAT (4F3.1) 

READ 500 (Z(M),M=1,4) 

DO 16 1=1,3 
DO 15 M=1,4 
DO 14 J = 1,10 
DO 13 K=1, 17 
CM1N1=3. 

CMAX1=1.. 

DO 12 L = 1,4 

A = (R1 ( J )»R2(K)+2. 1 1E3*Z(M )»R2(K)*E XPF ( .0183*T( L) ) )/ 1 .E6 
E = 2. 1 1E3*Z( M) *EXPF( .01«3*T(L ) )»(V(I)+(2.*R1(J))/1.E3)+V(I)»R1(J) 

F = 8.904E3*Z(M) *Z( M) *EXPF( .0366*T(L ) )+R2(K) *( V( I )-300. ) 

B = ( E + F )/ 1.E3 

C = 4.22*Z (M)»EXPF( .0183*T(L ) ) *( V( I ) -300. ) 

CTWO = I-B+SQRTF( (B*B)-(4.*A*C) ) )/(2.*A) 

IF (CTWO-1.5) 13, 13,8 
. 8 IF (CTW0-2.5)9, 13, 13 

9 CMAX1=MAX1F(CMAX1,CTW0) 

CMIN1=MIN1F(CMIN1,CTW0) 

12 CONTINUE 
DELTA=CMAX1-CMIN1 

600 F0RMAT(F5.0,2F13.0,F8. 1,F13.5) 

WRITE OUTPUT TAPE 4 , 600 , V ( I ) ,R 1 ( J ) ,R 2 ( k ) , Z ( M ) , OE LT A 

13 CONTINUE 

14 CONTINUE 

15 CONTINUE 

16 CONTINUE 
END FILE 4 
PAUSE 1 
END 

END 

60. 80.100. 

40000. 50000. 60000. 70000. 80000. 90000.100000.110000.120000.130000. 
200000.220000.240000.260000.280000.300000. 320000 . 340000. 360000 . 3« 0000 . 400000 . 
420000.440000.460000.480000.500000.520000. 

25.35.45.55. 

2.51.71.0 .4 
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PROGRAM COMPA 

DIMLNSION V( i) ,RJ,(,10) ,R?( 1 7 ) , T ( it ) , 7 ( 4 ) 

100 FORMAT! JF4.0 ) ‘ 

READ IOO(V{ I ), I=(l,3) 

200 FORMAT ( I0F7.0 ) 

READ 200 (Rl(J), J= 1 , 10) 

300 FORMAT ( 1 IF 7.0) ' ' 

READ 300(R2( K ) ,K= 1,.I7 ) 

400 FORMAT {4F3.0) 

READ 400 (T(L), L=jl .4.) 

500 FORMAT ( 4F3. 1 ) ’ 

READ 500 (Z(M),M~1,4) 

DO 16 1=1,3 
DO 15 M=l,4 
DO 14 J = 1,10 
DO 13 K= 1 , 1 7 
CMIM 1=3. 

CMAX1=1. 

DO 12 L = 1,4 

OCTWO={ {300.-V{I))/(Rl(J)+2.IIE3»7(M)«EXPF(.01H3*T(L)))-(V([)/ 

1R2(K ) ) )* I .E3 
IF {CTWO-1.5) 13, 13, n 

8 IF (CTW0-2.5)9, 13,13 

9 CMAX 1=MAX 1F( CMAX I , CTWO ) 

CMIN 1=MIN 1 F( CMIN I ,CTWO) 

12 CONTINUE 
DELTA=CMAX l-CMIN 1 

600 FORMAT (F 5.0, 2F I 3 .0 , F8 . 1 , F 1 3 . 5 ) 

WRITE OUTPUT TAPE 4 , 600 , V ( I ) , R 1 ( J ) , R 2 { K ) , Z ( M ) , DE L T A 

13 CONTINUE 

14 CONTINUE 

15 CONTINUE 

16 CONTINUE 
END FILE 4 
PAUSE 1 
END 

END 

60. 80.100. 

40000. 50000. 60000. 70000. 80000. 90000.100000.110000.120000.130000. 
200000.220000.240000.260000.280000.300000. 320000.340000.360000.380000.400000. 
420000.440000.460000.480000.500000.520000. 

2.51.71.0 .4 
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rair. ^.rogra^.- ill'iistratas one of the prto: .jr_’ .untages 

of the digital computer In evaluating circuits which have known para- 
meters of the type requiring incremental analysis. The flow chart of 
the main program called neon or light which uses the hyperbola or second 
order approximation is illustrated in figure 12, Here again every possi- 
ble combination of circuit parameters are systematically tested to deter- 
mine if they meet specif teat ions « With each set of design parameters 
the known parameters requiring increments are combined In every possible 
manner. If every combination of these incremental known parameters pass- 
es all logic tests with one set of design parameters, then this combina- 
tion of design parameters is recorded on magnetic tape along with the 
maximum and minimum neon current for the particular voltage. This signi- 
fies that the present design parameters being considered have produced a 
circuit which has successfully passed all of the logic tests for the un- 
predictable known parameters Z and T, The FORTRAN programs for neon and 
light are depicted In figures 13 and 14, 

Each set of output data in figure 15 is a successful design combina- 
tion. In other words for a neon tube maintaining voltage of 60 volts, 
a series resistance of 90 [< ohms, and a shunt resistance of 200 K ohms^ 
the neon current will only vary betv/een 2.32944 and 2.01816 milliamperes 
for the hyperbola approximation and between 2,32938 and 2,01802 mlllamper- 
es for the second order approximation. This variation in neon current in- 
cludes any variation in the characteristics of the photoconductor due to 
temperature or light intensity within the assumed possible limits. The 
last two columns of the output data, maximum and minimam neon currents for 
the chosen design parameters, have no immediate use in the synthesis of 
the circuit. Since these currents had to be calculated to evaluate the 
success of the circuit. It was an easy task to store them in memory and 
print out the maximum and inlnlirium values for each design parsiuieter set® 
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Fig. 12 
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PROGRAM NEON 

DIMENSION V( 3 1 ,R 1 ( 2 I 1 ,R2 ( 341 ( 5) 

FORMAT! 3F4.0 1 

READ 100(V( 1 ), 1=1,3) 

FORMAT ( 1 1F7.0) 

READ 200 (Rl( J),J=1,21 ) 

READ 200 (R2(K),K=1,34) 

FORMAT (5F3.1) 

READ 500 (Z(M),M=1,5) 

DO 16 1=1,3 
DO 15 J=1,2I 
DO 14 K=1,34 
CMIN 1 = 3. 

CMAX 1= I 



CMAX 1= 1 . 

DO 13 M=1,5 

A1 = (R1(J)*R2(K)+2.11E3*Z(M)»R2(K)*EXPF( .4575 ) 1 / 1 .E6 
E 1=2. 1 1E3*Z ( M)*EXPF( . 4575 ) * { V( I) + ( 2. *R 1 ( J ) 1 / 1 . E3 1 + V ( I) *R I ( J ) 
F1=6.,904E3»Z(M )*Z(M)*EXPF( . 9 1 50 ) +R2 ( K ) * ( V( I 1-300. 1 
B1=(E1+F1)/1.E3 

C1=4.22*Z(M1*EXPF(.4575)*(V(I 1-300.1 
CT1=(-B1+SQRTF( (B1*B11-(4.*A1*C)111/(2.*A1) 

IF(CT1-1.51 14,4,4 
IF(CT1-2.5) 5,5,14 

A2= ( R 1 ( J 1*R2(K W2. 1 1E3*Z(M)*R2(K)*EXPF ( .6405 1 1/1 .E6 
E2 = 2. 11E3*Z(M1*EXPF( .6405) »(V( I 1 + ( 2. *R 1 ( J 1 1 / 1 . E3 1 +V ( I 1 *R 1 ( J 1 
F2 = 8.904E3*Z(M)*Z(M1*EXPF( 1.281) + R2(K1»(V(I 1-300. 1 
B2=( E2+F21/1.E3 

C2 = 4.22«Z( Ml *EXPF( .64051*( V( I 1-300. 1 
CT2=(-B2+SQRTF( ( B2*B2 1- ( 4 . *A2*C2 1 11/(2. 

IF(CT2-1.51 14,6,6 



»R1( J) 



CT2=(-B2+SQRTF( ( B2*B2 1- ( 4 . *A2*C2 1 1 )/(2.*A2) 

IF(CT2-1.51 14,6,6 

6 IF(CT2-2.5) 7,7,14 

7 A3=( R1 ( J )*R2(K W2. 1 1 E3*Z ( M ) *R2 ( K 1 *EX PF ( . 8235 1 1 / 1 . E6 
E3 = 2. 1 1E3*Z ( M)»EXPF( .8235 1* ( V( I 1 + ( 2J *R 1 ( J 1 1 / 1 . 63 1 + V ( 1 1 *F 
F3 = 8.904E3*Z(M) *Z( M)*EXPF( 1 . 647 1 +R2( K 1 * ( V ( I 1-300. 1 
B3=( E3+F3)/ 1.E3 

C3 = 4.22»Z(M1 *6XPF( . 8235 1 * ( V ( 1 1 -300. 1 
CT3=(-B3+SQRTF( ( B3* B3 1 - ( 4 . *A3*C3 1 1 1/(2.* A3) 

IF(CT3-1.51 14,8,8 
IF(CT3-2.5) 9,9, 14 

A4=(R1( Jl*R2(K) + 2. 1 1 E 3*Z ( M 1 *R2 ( K 1 * EX PF ( 1 .0065 1 1 / 1 . E6 
E4 = 2. 1 1E3*Z(M)*EXPF( 1.00651*(V(n+(2.*R)(Jll/1.E31+V(Il*Rl(J) 
F4 = 8.904E3»Z(M1 *Z(M1»EXPF(2.013)+R2(K1* ( V( I 1-300.1 
B4=( E4+F41/ 1 .E3 

C4 = 4.22*Z(M)*EXPF( 1.0065 1*(V( I 1-300. 1 
CT4=(-B4+S0RTF( ( B4*04 l-( 4.*A4«C4) 1 1/ (2. *A4 1 
IF(CT4-1.5) 14,10,10 

10 IF(CT4-2.5) 1 1, 1 1,14 

11 CMAA=MAX1F(CT1,CT2,CT3,CT41 
CMIN=MIN1F(CT1 ,CT2,CT3,CT4) 

DELTA=CMAA-CMIN 
IF(DELTA-. 15 1 12,12,14 

rMAVl— lUlAVICrfrMAA rMAVl ' 



8 

9 



IF(DELTA-. 15 1 12,12,14 

CMAX 1=MAX1F(CMAA,CMAX1 1 
CMIN 1=MIN1F( CMIN,CMIN 1 1 
CONTINUE 



CONTINUE 

FORMAT! F6.0, 2F 13.0, 2F 13. 51 

WRITE OUTPUT TAPE 4,600, V ! I'l^R 1 ! J 1 , R2 ! K 1 , CMAX 1 , CMI N 1 
CONTINUE 



CONTINUE 
CONTINUE 
CONTINUE 
END FILE 4 
PAUSE 1 
END 
END 

60. 80.100. 

40000. 50000. 60000. 70000. 75000^ BOOOO. 82000. 84000. 86000. 88000. 
92000. 94000. 96000. 98000 . 1 00000-. 1 05000. 1 1 0000 . 1 1 5000 . 1 20000 . 1 30000 . 
40000. 60000. 80000. 100000. 120000. 140000. 160000. 180000.200000.220000. 
260000.280000.300000.320000. 340000. 360000. 380000.400000.420000.440000. 
480000.500000.520000.540000.560000'. 580000. 600000.620000.640000.660000. 
700000. 

.4 .91.41.92.5 



Fig. 13 
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niMfcNS ION V( 5 ) ,K U P 1 ) , ( 7i4 ) ,Z( 5 ) 

FORMAI (51-4.0) 

RhAlJ 100 ( V( 1 ) , 1=1,5) 

format ( 1 IF/.O) 

RFAD POO ( R 1 ( J ) , J= 1 ,2 1 ) 

REAU 200 (R2(K) ,K=1 ,i4 1 
FORMAT (5F5.1) 

RFAD 500 (Z(M),M=) ,5) 

DO 16 1=1,5 
DO 15 J=1,2I 
DO 1 4 K= 1 , 34 
CM IN 1 = 3. 

CMAX 1=1. 

DO 13 M=l,5 

(JAl=?.ll*Z(M)*EXPF(.45^5)»(1.42ttR2(K)+5T0.P«/(M)»bXPF(.45Z5)f.P5118 
1*R1(J))+1.42*R1(J)*R2(K)/I.t S 

EI=2795.6*Z(M)*EXPF(.45Z5)»(Rl(J)/l.E3+2.1l«2(M)«EXPF(.4575)-(300. 
1-V ( 1 ) )/5.273)-R I ( J ) »R2( K ) / 1 ,t 5 

Fl=2. 11*Z(M)«EXPF( .4575)*( 1.4?E3»V(I )-R2(K) ) + l.42*(V( I)*Rl(J)--500. 
I*R2(K)+V( I )»R2(K) ) 

R1='E 1+F 1 

C 1=Z ( M) »EXPF( .45 75 ) *( P 7V5.6<^ ( V ( I ) - 300 . ) -P . 1 1 h3< 

1V( I ) ) + (300.-V( I))*R2(K)-V(I)'^RI(J) 

CT 1=(-H USORTFI ( 111 *R 1 ) - ( 4 . * A 1 1 ) ) ) / ( 2 . * A 1 ) 

IF(CT1-1.5) 14,4,4 

IF(CTI-2.5) 5,5,14 

A2 = 2.1l*Z(M)»bXPF( .6405 )*( I . 42 »R2 ( K ) +53 0 . 2 « / ( M ) « b XP T ( . 64 05 ) + . 25 1 1 8 
1*R 1 ( J ) ) t 1 .42»R 1 ( J )*R?(K ) / 1 .-F 3 

E2 = 2 795.6*Z (M)*bXPb( .6405) (R 1 ( J) / 1.F3 + 2. 11 «Z(M) *£XPF (. 64 05 )-( 300. 
I-V( I ) )/S.273)-Rl(J)»R2(K)/l .E3 

F2 = 2. 11»Z(M)*EXPF(.6405)*( I. 4PE 3»V ( I)-R2(K)) + 1.4P»(V(l)»Rl(J )-300. 
1 «R2( K) +V( I ) »R2(K ) ) 

R2=E2+F2 
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1 1 
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13 

600 
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16 



C2=Z ( M) *EXPF( .6405 ) *( ?795.6« ( V( I )-300. )-2. 1 1E3* 

1V( I ) ) + ( 300.- V( I ) )*R2( K )-V( I ) *R 1 ( J) 

CT2=(-H2+SQRTF( ( BP*B2 )-( 4 .* A2»C2 )))/(?. *A2 ) 

IF(CT2-1.5) 14,6,6 

6 IF(CT2-2.5) 7,7,14 

7 A3=2. ll»Z(M)»EXPF(.B?3S)*(l.4P*R2(K)+530.2»Z(M)»EXPF(.8235)t.25ll8 
1*R1( J) ) + 1.42»R1( J)*R2(K)/l.b!3 

E3=P795.6«Z ( M ) *bXPF ( .023 5 ) » ( R J ( J ) / 1 . E3 + 2 . 1 1 * Z ( M ) *EXPF ( . 82 35 ) - ( 300. 
I-V( I ) )/5.273)-Rl ( J) »R2(K)/,I.,L3 

F3 = P. 1 l*Z( M) •bXPb( . 8235 ) »(M .'4P,E3*V( i )-R2 (K) ) + 1 .42*( V( 1 ) »Rl ( J )-300. 
1*R2(K)+V( I )*R2(K) ) ‘ ' 

B3=b3+b3 

C3 = Z ( M)*EXPr ( .8235 )*( 279, 5.6« ( V ( I ) - 300. ) -2 . I 1 F3* 

1V( I ) ) + ( 300.-V1 I, ) )*R2( K )-VU ) «R 1 ( J) 

Cr3=(-R3 + SQRTP( ( R3*B3 )-('4.. *a3*C3 ) ) ) / ( P. *A3) 

IF(CT3-1.5) 14, 8,0 

8 IF(CT3-2.5) 9,9, 14 

9 A4 = 2. 1 1*Z( M )*£XPF( 1 .0065i)*,( 1 .42*R2(K) +530. 2 » Z ( M ) *E X PI ( 1 .0065 ) *-.251 Id 
1*R1( J) )+l.42*Rl( J)*RP(K)/1.I 3 “ 

E4 = 2795.6»Z ( M ) »tXPF ( 1 . 0065>) »'(R H J ) / 1 .£ 3 + 2. 1 1 »Z ( M ) *£XPF ( 1 .0065) - ( 30 
I 0.-V(I))/5.273)-RI(J)*R2(K)/l.b3 

b4 = 2. 1 1 «Z( M ) «EXPF( 1 .006 5 ) Hi . 42E 3 « V ( 1 ) -R2 ( K ) ) + 1 . 4 2» ( V ( I ) » 

IR I ( J )-300. *R2( K ) + V(<n*R2 (K1 ) 
rU=F4+F4 ' ' 

C4 = Z (M) »EXPF( 1 .0065 )« (2795.6*( V( I )-300. )-2. 1 ll:3« 

1V( I ) ) + (300.-V( I ) )«R2(K)-V< I )»RJ(J) 

CT4= ( -R4 + SQRTF ( ( H4*B4 ) - ( 4. *A4».C4 ) ) ) / ( 2 . *A4 ) 

Ib(Cr4-).5) 14,10,10'' 

IF(CT4-2.5) 1 1, 1 1, E4 

CMAA = MAX IFICT l,CTP,CT3,CT4 ) 

CMIN = MIN1F(CT 1 ,CT2,CI3,C 14 ) 

DEL r A=CMAA-CM1N 
IF (DELTA-. 15) 12,12,14 

CMAX 1 = MAX 1 F ( CMAA,CMAX‘l 1 
CMIN1 = MIN1F(CMIN,CMIN1 ) 

CONT INUE 

FORMAT! F4.0, 2F 1 3.0, Pr 13.5) 

WRllE OUTPUT TAPE 4,600, V ( 1 ) , R 1 ( J ) , R? ( K ) , C M AX 1 , CM I N 1 
CONT INUb 



CONTINUE. 
CONTINUE 
END FILE 
PAUSE I 
END 
END 



Fig. 14 




I 




VfA. 


Rt. 


Re. 


Ckax i 


Cmin 1 




(»0 . 






1 . M ) f 


1 . / 




/)0 . 






1 . \ f 


1 . ( 




6l) • 




1 ' ) . 


'' . '' / )l 1 . 


1 . '4/8 




60 . 


1 ' ' , , 


18; 


. ) ' ] > 


1 . / / '6- 




6 0 . 


' i . 


P. '' 6.. ' 


. ', \ yf 


1 . / 8 ) K 




6 0 . 




1 (icr ! ( . 


1 )*6'i 


1 . 4 / / 8 !) 




(f 0. 




I ■)>' , 


P ^ .8 1 


8 . ;, 8t 0 




60, 




/ ' ' ( ' 0 0 . 


’ . > ^ 


8 . ' /,8 6 ; 




60. 






’ . I •; 1 ; 


6 ’*8 




I'l u . 


, , 


, .')( . 


.8 807 


. 1 76 , 1 




60. 




;0 6. 


. )| ')■' 8 8 


’ . 1 > ) 8 j 




60. 
60 . 


. . 

O'" ' : . 


;0i-. 


^ . -1 0 ■; 1 0 
/.It n( 


8 . U 8 
. 18 1)8 


R?OG(?AfA 


oO . 
60. 


. 

^ 0 ) . 


^ >;;'!()!• . 
6!)0;)( . 


/ . '1 ■/ /() ■ . 
1 . 0 ■ , ^, 7 0 


. 1 6 P) 1 
1 . 8 8r,8 8 


MEo/nI 


60. 


'/{){ )l' . 




1 . ; /( M 


1 . 6>7)o / 




60. 


. 


1 ' a* - ' ) f , 


6 . 6 ') '))♦ ) 


1 . 8 M 8 / 




60. 


•/■or M, , 


1 r c' ) 1 ,1 . 


8.1 / ^ ; 


1 . 8 8 8 / 




6 0 . 


V(M:0 . 


1 ■.v;000 . 


6 . . >'6 ') 


1 . / 11 'U; 




60. 


/Oil ;ij . 


100(*0- . 


8 . v-,8/^ 


1 . >84 V 




60. 


.^0 ' 0 . , 


1 '.oooc. 


8 . ■' / 6 8 


1./ 867 




60 . 


/0(iUc . 


'•'■.'Cm): . 


7 . , )4 M 


8 . 7 1 1 0 




60. 




V,Om)()(> . 


0 


8 . : 7, 1 r . 




60. 


/{■OOM . 


?U()()0(; . 


8 . ) 8 ' ^ ; 


^ . 76 1 / / 




6 0 . 




.'JO . 


8. :7 8 


8 . 0 8/08 




60. 


7! r.iWi. 


. 


; .'1 1 1(77 


/ • / 1 




60. 


VGwft'/ . 


■)0'^('0 0 . 


7 . I< / ■ 0 6 


8 . 106> 8 




60. 


vonoa . 


•P' f*Of)() , 


7. 'Cl O' f 


. 1 188 / 


’ 


60 . 


7 COO' . 


0( 0 0 . 


.4'. 1 


. 1 ' ii 8 






'60. 


. 


6 11 00 t,: . 


1.81 8 Dm 


1 . 87(^78 


" 


60. 


86,0 0:') . 


0 00 0 . 


1 . 7^4 9 18 


1 . 8 8 8 8 1 




60. 




1 00000. 


8.0 >87 8 


1 .68088 




60. 


680 0\) . 


1 /’OOOu . 


8 . 1 > 1 )< S 


1 . ‘^7 >6 8 




60 . 


H n 0 ^ . 


1 1| 1 00 0 . 


8.8') 1 '1 8) 


1 . 7 58 8)) 




60. 


8r0d:) . 


IGOOOO. 


7 . 8 1 )M 8 


1 . 7 96, >0 




60. 


86 0 00 . 


1 soooo. 


8 . 8 8 8 8 8' 


8.0 ^86^^ 




60. 


8^ 00 J . 


7 0000 0 . 


8,8 6 


8 . o6 8i*u 




60. 


rU-(U)'j . 


7 7(1 ‘10 0. 


8 . U 1 1 


8 . 0 66 1 9 




60. 


'OUUm) . 


71|C.'UJ 0. 


8 . U ',768 


8 . 106 




60. 


H ) 0: '' . 


70000C). 


8 .’*86 81 


8 • 1 8 850 




60. • 


■•’800 . 


7. -1 0 00c . 


8.)‘ 8 8^6 


8. 1 py 




60. 


86000 . 


OfU.OOO . 


8 . )»678)| 


8 , 1 5 1 )i (j 




60 . 


68000 . 


<70000 . 


8 . )J >766 


8 . 16 866 




60. 

60. 

60. 


vOO 0 0, 
8000' . 
VOOOC . 


6 0000. 
■•'.O'^Oo. 

1 .)00()!j . 


1 . 6 8 8 8 ! 
1 . ^’ > • ) S 
8 . 8 84 1 


1.88817 
I . 6‘'‘684 
1 . 8 4 0 1 


R^OG(?A^^ 


60. 


VOOOO . 


1 '■000 0. 


8 . 1 8 8 ) 


1 . 6 3 8 


LlGfAT 


60. 


8f)nofj . 


Hi 000 0. 


8. 8('888 


1.711 81 


60. 


900f''0 . 


lOliOOl). 


8. 8 /.DO 


1 . 7 84>i , 


60 . 


VO 0 00 • 


1 i 0i)0 0 . 


8 . 8'v’6*{6 


1.76 7 5 <\ 




60. 


V('0 00 . 


700000 . 


8 . 8 > > 8 


8.0 1 'D)/ 




60. 


voono , 


770000. 


8 . ' 6 •; 8 


8 . OJ) 1 8.' 




60. 


VOOO'i . 


7'ii:'ioo. 


8. (> 78 


. M/, ] 8(, 




60 . 


/I'OO.'' . 


7oiH)U0 . 


8 . 8 M 


P n 8 -0 1 




60 . 


voono . 


7 O'lOO. 


8, *M )|7 8 


8 . -)7 3 8r, 




60. 


y()(M) ; . 


■I 0(1 0-0 0 . 


8 . U 86)i 8 


’ . 1 ( 7 8 .3 




60. 


70000 . 


'0000 . 


8. )J46 >’i 


’.11 18 




60. 


;of)0o . 


^’l(.UOC . 


8 . 4') 1 0.: 


8.1/' 38 




60. 


voooo . 


v.iii'O 1 , 


8. 1 


8 . 1^8 * / 




60. 


8000'* . 


■1 ;oouv). 


P ’1 8 0 >6 


. 14 V >4 




60. 


;ooot. . 


ii I'oOl' (■ . 


r " 18 


^ . 1 7 >1 




60. 


• /'('OO') . 


li 1 . 


> . '1 ■ M . ■ 8 


^.18 >-( 




6(J . 


;^C0 Of* . 


h 'i r ( t'l '■ , 


'.•i > 1<- : 


>. i ^88( 




60. 


6 (0)0 . 


ii />( 1 Ol) . 


■' . 4 .M" ’ 


' . 1 8 1 / r 




60. 


; / . u) ' 1 . 


■,l 7> 


1 , • 8 1 


1 . 1 ' ' 





4 0 



Fig. 15 
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These iuniukani neon €n,rrents -'h:: have value tox the 

after the design has been choi^an^ and all that tBiue\±vc is the datermtna* 
tion of the particular r/ialntaining voltage in any one circuit* This will 
be of immediate us<e provided that the neon tube voltage is close to one 
of the increments of voltage used In the digital computer program; in 
other cases, interpolation is req?jired« 

The digital computer program which Is flow charted in figure 12 and 
represented in FORTR/IN statements in figure 13 has a computer running time 
of 8 minutes and 30 seconds* This output data still requires interpret^® 
tion, evaluation and presentation in a useable manner before it can be 
considered as a trtie design asset* 
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xUjx t€ e'^tei.-jiivt^ irtevprecst ion of the •■ii^ital eosiiptuter 
printo.U', it Is i 4^^1t a ivantageoys tu considst the n^on photoconductor 
circuit OP slrj let ^asls* l..^t un 'otake the follootitg assu:*^pt ions : 1) 

The Incoming VDlt«^ge of the neon tube ts fixed at 60 volts* 2) The 
knoWi paraMCter has only one volue 2*5* 3) Temperature varia- 

tions are such t\>\t the amhient tai'.paratur*i will vary between 25‘^ C and 
55"^ C, vvlth these three assumptions tv;o circuits will be considered, one 
wuth ccTiipensatlon and on-e without tha conipensating shunt photoconductor* 
The printout for the compensated and m^compensated circuits are sho\^ in 
figure 16, These two printouts are of the same format as the main program 
except that the maxin'.um and alninum currents are now solely attributed to 
temperature variations* 

A plot of the de:ji 3 n parameters R2 versus Rl^ which meat output 
specif lest namely rieicn curr^^nt between 1,5 and 2,5 mil llaraper^^s, 

is constructed for both the co^Kpensated and uncofi.pcnsated photoconductor 
circuits, Th-ese plots ar» figures 17 and 18 respectively. The uncom- 
pensated photocondactor circuit has 5t minimum value of shunt resistance* 
This constraint is brought about by the ph^^slcal requirement of a suf- 
ficient voltage drop across tha shunt resistance to fire the neon tube* 

The compensEted photoconductor circuit has a dark photocondiu^tive re* 
sistance In series with the shunt resistance prior to firing which jLs in 
chc u.c 2 :ioh » region. Therefore, wn^y the uncompensated ::lrcui. requires 
a uirilmuxi ^alue ol shu?.t leslstancc. ..ssuiLiag a iicireg voltage of 135 
olts which is sn ohsev^red aveisga value, Jii 2 u the follov.’^vig ca ici<ila- 
ti ‘US for the minimum value of R2 at L 2,5 apply t 
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effect of the r uiient due only to tei pcrature variations* 

These nun-bets vithlr' th.- ate'^ represent tlie difference between the maximuru 
and minimum neon cvarent in mil liamperes for that particular set of design 
{parameters (R1 and R2), A comparison of the two plots reveals that compen- 
sation is actually achieved by insertion of the shunt photcconductor. A 
second conclusion that can readily be made is that for a particular value of 
R1 in the uncompensated area plot the variation in current is constant and 
independent of the value of the shunt resistance R2* This constant value 
is also the limiting value v^^hich the compensated area plot approaches as 
the shunt resistance becomes Infinite. Analytically this last result 
could have been foreseen since for large values cf shunt resistance the 
effect of the compensating photoconductor on the total series resistance 
diminishes and eventually approaches zero. This simple neon photoconductor 
circuit will be employed later in this section as an illustrative example 
of best design procedure. 

The digital computer printout of the main program whether it is the 
hyperbolic or ’^second approximation is of the same format (See 

Figure 15 for printout format). For each value of series resistance Rl, 
those values of shunt raslstarsce F2 which it^ark the upper and lower limits 
of successful circuit combinations in the printout are plotted on graph 
paper producing areas. In the case of the hyperbola only two of these 
plots are necessary, one for the lowest maintaining voltage encountered 

V ^ 

and one for the highest ^maintaining voltage since — - is always nega* 

tive. The **second order’" approximation expression is nore In'volved and 

it cannot bs mathematically proven to have a negative for all 

V 

design sets to be considered. In this case incrcinents of maintaining 
voltage v;erfe used Ir the iigltal ronputer program and as many plots as 
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u • , .-h t -v plot Is for 

a [•articui^s.. u^xiiita ir ^'.--.•.‘C3rj;vi, al’ Z’s and all 'tVs^ In both 

apx^roxiiTfat ions the resuiting vhether two at In the case of the 

hyperbola or many as in the secern i ardei, can now be overlapped to pro- 
duce a final area which is cotnmon to all of the individual areas. This 
final area If it exists has satisfied all of the restrictions i^fi|>osed on 
temperature^ 2, maintaining voltage and output neon current. The resultant 
area for the hyperbolic approximation is shown in figure 20. 

In order to observe the effects of irore rescrictive specif ications, 
the printout data of both the hyperbola and ^second order*’ apj roxiniations 
was plotted using those upper and lower limits of the shunt resistance 
that yielded an output r%eon current betiveen 1.8 and 2,2 mllliamperes. 

After having plotted the Individual area plots^ the overlapping pruceduie 
was employed to find tha coiV4fi;on area. A comracn area Is non^exlstant . If 
this specification of output currant had been used vice 1.5 - 2.5 iritni- 
amperes, then the absence of an area indicates that there is no design para» 
meter combination that Is satisfactory for the entire maintaining voltage 
range. 

The amount of compensation effected can ba seen by analysis of the 
difference in maximum and minimum neon currants due only to temperature 
variations. Refer to the simple example ment letted previously where the 
maintaining voltage was fixed at 60 %^olts and parameter Z fixed at 2.5 
(See figure 17). The mapped points in the interior of the area plot in 
figure 17 represent the difference between the maximum and minimum neon 
currents in milliamperes due solely to temperature variation for that 
particular set of design parameters (FI and R2), Inspection of the plot 
reveals that the deviation due to rertiperat are for the same value of R2, 
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■ " a r 35? ^ t r »■ i < T 

stin • th** *.r,|j,-.-i ^ sair-e v^^.ue of ^ the 

devlatlcti) .aa to ter- f eiatut^ h^s its srr^t^llest value vhf:'i! R2 tahes on its 
smallest j.ii.sslb1e valoe fTisi.le the plotted are&* The smsie conclusions 
are true for other similar plots. To se^ the overall effect on the riain 
program for the hyperbola and ^^seccnd order*' approximations, the dlffcirence 
in neon currents due to temperature variations (called delta in the compater 
program) was restricted to first *015 and then ,010 millla^riperes^ The re- 
sultant area plots for the hypertola without delta and with delta equal to 
,015 are shown in figure 21. With the program rerun using delta of *010, 
a resultant area did not exist. A comparison of the final area curves with 
and v/ithout delta reveals that the left edge of the resultant areas which 
employed the delta restriction are moved further to the right, where as 
the right edge is unaltered. 

The resultant area curves^ since they are in effect the common area 
of the individual areas, cannot be mapped for the variation in neon currents. 
Any point in the resultant area meets all of the specifications that were 
required and is a good design set. To determine the best part of the^ final 
design area we will resort to inductive reasoning. In inductive reasoning, 
a set of individual cases is studied by the experimental method, and from 
the observations made, a general principle is formed. The principle formed 
by Inductive reasoning is reliable only when all possible instances have 
been examined. In the Individual areas, the mapping procedure resulted in 
the conclusions that the best design set would be the smallest value of R2 
and the largest value of R1 that fell within the plotted area. Applying in- 
ductive reasoning to this information, the best set in the resultant area 
curve is also the smallest value of P2 and the largest value of Rl. Consider 
at ion also has to be given to the variation in the design parameters Rl and 
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R2 , * j *,rt then the, afex cE the ^uea would be the best 

design set. R1 R2 uiy ha/e tolerar.ces themselves of from 5 to 10 per 
cent of theit" oonupal value, hence some point that could be the center of 
a circle whose radius Is the resistance tolerance and where the circle fell 
entirely within the resultant area at the lowest part, would be the best 
design set. 

The resultant area curves of the hyperbola and second order approxitna^ 
tions are almost identical in position and shape. Minor variations in shape 
can be attributed to the plotting of the individual computer incremental out- 
put data. As an example of this problem in plotting consider a design in- 
crement change of 2 kohms in R1 where the value of the design parameter re- 
sulted in a successful circuit prior to the increment change, but failed 
after the increTnant change. Theoretically thon the plotted curve could pass 
through the first point and it also could pass through any point to the right 
of this original point up to but not including the value of the original 
point plus the increment. Here then is another factor in considering the 
number of design sets to be investigated versus computer time. This factor 
should not carry too much weight since in the final design, the designer 
will pick a point somewhat removed from the actual boundary curve itself. 

The resultant areas as obtained by the two approximations should be the same 
for large values of R2, The difference is noticeable v/here R2 is 280 kohms 
or less. In this range of R2 compensation is more effective and the second 
order area is a more accurate representation. However, In the final analysis 
the differences between the two approximations prove to be very minor. 

Returning to figure 1*^ where the voltage is fixed at 80 volts and Z is 
2.5, it appears that R2 at R1 s 83 kilohms. The digital computer 
printout data when plotted reveals that the lower curve is In effect the 
locus of combinations of R1 and P2 which result ix\ neon currents of 1.5 
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■ ^ ^ i ilai locus ufjt: Cull errs of 

mi llan-pcres* c jir> ctaa^. ly only r^pproxirate locli of these 

constant current \c\lucs since the joints which were used to construct the 

curves were obtalneJ by using inc^cr^^ntal values of design parameters* 

0 

The system e<iuation which was of the fonn AI‘ f BI r C - 0 \^as 
solved explicitly on the digital computer for as a function of five 
parameteres* If the system equation is rewritten so as to express R2 as 
a function of the remaining parameteres, the result using the hyperbola 
approximation would 



Using a neon current of 2*5 mil liairperes, a voltage of 60 volts, a 
Z of 2,5 and a iiiean temperature, of 40 the value of R1 to make the de- 
nominator of the equation 2^ero is 35 kohms* Using this equation similar 
results could be predicted for other maintaining voltages, namely that in 
the plot of design paramter R2 versus design parameter HI there arc values 
of R1 where R2 goes to Infinity. Since the resultant area is constructed 
from the individual area plots we cannot expect to have a closed resultant 
design area plot. 

The insight gained from the interpretation of the digital computer data 
for the circuit synthesis problem revealed that the constructed boundary line 
are in essence locii of constant neon currents, namely the extremes of the 
specified neon current limits of 1*5 anc! 2.5 millismperes. This suggests an 
alternate technique for solving the circuit synthesis problem. This alter- 
nate technique lends itself to solution using an analogue computer. 
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7* u -ic toMc (.omputer Technique 

The technique t o be loved with the analogue corputer parallels the 
development of the digital computer technique to a certain extent. Known 
and design parameters are as defined previously. The resultant circuit 
equation is manipulated to express cne design parameter explicitly as a 
function of the remaining design and knovm parameters* Increments of the 
known parameters and end points of the known parameters are determined as 
in the digital computer technique* The explicit circuit equation is then 
converted into an electronic nnalcgue* The electronic analogue computer is 
set up so that its equations have the same mathematical form as the equation 
for the physical system under considerat ion* In the computer, voltages are 
used to represent various physical quantities such as velocity, position, etc. 
Scaling factors are used to convert the physical quantities into appropriate 
numerical values in the computer* The mathematical functions of addition, 
multiplication, etc., are performed by electronic operational amplifiers 
using scaled voltages as inputs. 

large variations in knovTi pi?ran<eters can be readily undertaken pro^’^id- 
ed the range in vnri^tlor is considered during the scaling process* Incre- 
mental charges in the knoiym parameters are accomplished by merely varying 
the coefficient potent io^i eters that are associated with the parameter urder 
Investigation. Badical changes in the known parameters may have to be 
handled by varying the components of the operational amplifiers, or the 
insertion of one or more additonal amplifiers. 

A design parameter that is implicit in the circuit equation is varied 
by the application of a raiMp voltage input causing the explicit parameter 
output voltage to be the continuous solution for a fixed set of known para- 
meters. Utilizaticn of an X-V plotter will present a graphical representation 
of a design paranaeter voltage as a function of another design parameter voltage. 
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tHfi X** X plottiir* ftjiiiT'bjr u n p-j^r-fi^T'e^twcs thnt. cm'* b© vstisd on sny 

one sheet o£ the riotAer becc'wjs avlJT'.et Vahen o less of identity In the 

plotted curves results. This pro t>4ure is repented until graphical plots 
are obtained lor every ineremert of each of the known paraitseters. These plots 
are antiscaled on the X-Y plotter so thst the graphical plots which are ob- 
tained in this manner are design parameters as used in the original physical 
system prior to the analogue conversion. Each plot has only investigated one 
or more variations in known pararoeters. If all of the plots are overlayed, 
the area x4Uch is comiron to the acceptable area of all the graphs is the final 
design area desired. Any spoint In this final coranon area meets all of the re- 
quired specifications. 

Rearrangement of the basic circuit equation employing the hyperbola 
approximation for the pbotoconductors resulted in the form: 
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The neon current I is now an end point type of known parameter rather 
than the output quantity as in the digital computer ter.hnique. The two neon 



current values of interest are 1.5 and ?,5 milliamperas. Physical quanti- 
ties such as ohms. ampereSj, and maintaining voltage in the circuit equation 
are transformed intc voltages for use in the electronic analogue by appropri- 
ate scaling. Known parameters of temperature and Z are divided intc in- 
crements as before. A change in the value of any known parameter is accomp- 
lished by altering the individual potentiometers that are associated with 
the parameter under investigation. Since the range of all the know para- 
meters vary by less than ? factor of ten» it is possible to utillta one analo- 
gue computer circuit aid accourt for all of the variations by suitable scal- 
ing of the physical parameters, r schematic of the analogue computer 
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ampllClci 'u zi^ute 22, 

A step voltag.^ npplio] to 1 r.»:er;^raL Ing amplifier F produces a 
ramp voltage ^'^hich corresponds to des igt^ parameter Rl, This continuous vari* 
ation of design paranetei: Rl, rith a fixed set of temperature, Z, maintain* 
ing voltage and neon current will produce a continuous solution for the ex-* 
plicit design parameter R2, Potentiometers a4, a5, a7, alC, all and al2 are 
held constant at one value throughout the problem. These constant potentio- 
meter settings are determined by scaling physical quantities which are in- 
variant. The setting of the variable potentiometers at any one instance 
depend upon the known parameter that is being investigated. Table #2 lists 
the depender.ee of each known parameter on its associated potentiometers. 



Table ^^2 



Variation 
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Potentiometers 


a3. 


a8j a 15, a 16 
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Potentiometer 
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Potentiometers 


a9. 


al3, al4 



Figure 23 

Each of the known parameters one at a time is systematically varied by 
changing their associated potentiometers. The tv;o voltages in the electron- 
ic analogue which represent design parameters R1 and R2 are automatically 
plotted by inserting these continually varying voltages as Inputs to an 
X-Y plotter during any one solution which corresponds to a particular knov/n 
set of parameters. The output of the X-Y plotter is a graphical representa- 
tion of design parameter R2 versus design parameter Rl. A typical output 
curve for Z (pbotoconductor resistance magnitude multiplying factor) corres- 
ponding to .4 appears as in figure 2'*, This curve includes the effect of 
temperature, maintainiTig voltage and neon currents. Solid lines represent 
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maintainiug voltnges of volts ri i dished lines represent (^aintnlning 
voltages of ICO volts* The three temperatures investigated are 25^C, 40*^0, 
and 55*^C which are t*ecorled from right to left in groups of three lines* 

For each maintoining voltage* the lover gn up corresponds to 1.5 mllliamperes 
and the remaining group correspomls to 2*5 milliamperes. Similar curves for 
other values of Z appear in figures 25, 2*5 and 27. 

As wavS mentioned previously in section four, the hnoim parameters, 
temperature (T), and pbotoconductor resistance magnitude multiplying fac- 
tor (Z), were established as incremental types of parameters since --- - — 



and 



i: 

j A 



could not be proven to be of constant sign for the entire range 



of variation. Zxamination cf figures 24 through 27 discloses that for the 
maintaining voltage equal to 60 volts, the three temperature curves of t!ic 
low’er group (neon current equal to 1.5 milliamperes) cross over each other. 



This cross over is indicative of a point of inflection in or stated 

J T 

In other terns the value of ^ changed sign at the point of temperature 

J "a 

cross over, By overlaying figures 24 and 25 a similar effect Xvould be noticed 

for Z, namely chat for one value of temiperature, maintaining voltage and neon 

current the overlayed curves correspording to two cUffex^ent values of Z 

would cross over. This cross cvex indicates where a point of inflertion ex= 

*\ r 

Ists for , Aai overlapping technique as applied to the digital coiniuter 

} ^ 

areas is used in finding the area coirimon to the individual areas, '''he re- 
sultant design area is figure 28. 'lince all of the points of inflection for 
temperature and pbotoconductor resistance iriagnltude multiplying factor lie 
outsi^’e the resultant design area, these quantities are in effect end nolrt 
type of Vnown parameters if the term as used hare applies only to the result- 
ant design area. 

The design area as obtained from tha analogue computer in conjunction 
with the X-Y plotter meets the requiremevits of circuit perfcrr.ancc with any 
set of known para:-<eter5; . Tha Mnes forming the boundaries of the area are 
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ital cu' •' .%'e» . v.r^st j. ■■. tai.-eier cor-btoat ion f.-iuKL-t 

ed fi'Obi this dc??Ign ntca as wa-=? obtained iroim the digital 
there is no aJdit- n.al in fo 1113 1 ion available about the resultarst 
upon which the designer can base a choice* This perhaps is its 
as compared to the digital computer technique* 



ji- determir^“ 
er drtA^ since 
design area 
sole weakness 



64 



I 



1 




. R 2 IN V<lo(^a\s 




I 




P . • .iS U -nS . 

The i-rol le'n cf lu^v. o clicult to meet* sj er IficatioTiS ■ h le 

using f.vn remet* vs which beve Wicu^ rengf^s of vari^ition can be soU7ed using 
either of the coroputet techniques described* The designer has merely to 
select design paraniefcers from a fin^. 1 area v;ith the assurance that the 
circuit will raet the output specifications* 

Both the digital and the analogue computer techniques provide adequate 
design data* Each method has advantages: there is little difference in the 
amount of tine Involved from start of problem until solution* This time 
naturally depends upon the designer’s experience and familiarity v;lth a 
particular computer* The accuracy is comparable in either method; while 
the digital computer is more accurate for any individual calculation, the 
analogue has the advantage of providing a continuous solution rather than 
an incremental one. The analogue computer, when used with an X*Y plotter 
provides the area curves directly without the plotting necessary to reduce 
the digital data, Hox^?ever, the digital compiater can be easily progranmed 
to provide additional data useful In analysis of a specific, designed cir- 
cuit* There were no logic circuits employed in the analogue solution, 
therefore the significant test for delta I was not made in this solutioiu 

For the specific circuit used to demonstrate these techniques, it 
should be noted that the scheme for temperature comperisaticn did achieve 
the desired compensation as reference to figures 17 and 18 will show. The 
additional labor involved in employing the more complex second order approxi- 
mation is not warranted since there is no significant difference between the 
final design areas obtained with the two approximations* Figure 29 compares 
the results obtained with the two approximations using the digital computer 
and the hyperbolic approximat ion using the analogue computer* 
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Additioiial Investigation shujld be coniv'ictel to extend tbe study to 
include an c^nalysis of the circuit’s transient behavior, A method of in-- 
eluding logic circuitry in the analogue circuit would permit the additional 
test for delta I to be conducted, A circuit Involving the selection of 
three design parameters presents an interesting problem. In this case, 
instead of acceptable area graphs, the result would be a three dlrnensicn- 
al acceptable volume plot of the design parameters. 
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MAT {EMATICAL DERIVATION LEADINr: TO SECOND ORDER EXPRESSION FOP 
P TOTCCONDirCTOR KESISTANCE-NEON CURRENT RELATIONS'^IP 

To obtain a mathematical expression which approximates the photoconductor 
resistance variation V7lth neon current shoxim in curve B, the relationship 
is assumed to be of the form: 

/\ P,-.L < J I (■ k"j, I t-' 

The following points, selected from the curve were used to determine the 
unknown coefficients A, B, C and D. 



I in tna. 1.50 1.80 2,20 2.50 

R in Kohra 6.71 5.10 3.88 3.37 
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Substituting in (4a) and siotplifying. 



.7:6 X: c ■ ( I-* 

7.2 6 X>-r 



i V , 't* - -V X 



i- ‘■ X 



for (' I 



r> 6 3 X 



C L 



h‘i! I X /<- //.. J i I 



i r 



k'p 



£- /' ’ j i A / ' 

L '^s -i X ./ i j 



s K <■■ 



2 . 



^ .' \ ■ J 

i.H'H f 



To incorporate the effect of tenperature, rewrite in the following form; 
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